Scaling of Dynamical Decoupling for Spin Qubits 
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We investigate scaling of coherence time, T2, with the number of 7r-pulses, n^, in a singlet- 
triplet spin qubit using Carr-Purcell-Meiboom-Gill (CPMG) and concatenated dynamical decoupling 
(CDD) pulse sequences. For an even numbers of CPMG pulses, we find a power law, T2 oc (n^) 76 , 
with 7e = 0.72 ± 0.01, essentially independent of the envelope function used to extract T2. From this 
surprisingly robust value, a power-law model of the noise spectrum of the environment, S(u) ~ cj _/3 , 
yields f3 = j e /(l - 7e) = 2.6 ± 0.1. Model values for T 2 (n 7r ) using f3 = 2.6 for CPMG with both 
even and odd up to 32 and CDD orders 3 through 6 compare very well with experiment. 
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A variety of solid state systems are emerging as ef- 
fective platforms for studying decoherence and entangle- 
ment in controlled quantum systems [T]-[4] . Among them, 
quantum-dot-based spin qubits have recently achieved 
sufficient control and long coherence times [TJ [2] that new 
information about the noise environment of the qubit can 
be extracted, complementing related work in nitrogen- 
vacancy centers in diamond [4], superconducting qubits 
[3], trapped ions [5 , and neutral atoms [6]. 

Dynamical decoupling in the form of a sequence of n- 
pulses [TUTU] functions as a high-pass filter, thus pro- 
viding information about the spectral content of envi- 
ronmental noise [3-5] ITTHT4] . For spin qubits, the ef- 
fectiveness of various decoupling schemes at mitigating 
dephasing due to nuclear bath dynamics has been well 
studied theoretically [T5HT9] . Much less is known about 
mitigating the effects of charge noise, which couples to 
the qubit via gate dependent exchange interaction and 
through spatially varying Overhauser fields pQ. When 
the decoherence time, T 2 , is short compared to the energy 
relaxation time, Ti, — which is the case in this study — 
both the envelope of the coherence decay as well as the 
dependence of T2 on the number of 7r-pulses, n^, de- 
pend on the spectral density of the environment, S(uj). 
Knowledge of S(u) inferred from such measurements can 
in turn be used to design optimal decoupling sequences 
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In this Letter, we investigate scaling of T2 with 
the number of 7r-pulses for Carr-Purcell-Meiboom-Gill 
(CPMG) and concatenated dynamical decoupling (CDD) 
sequences in a GaAs two-electron singlet-triple qubit 
[Fig. 1(a)]. The coherence envelope is reasonably well 
described by the form exp( — (td/T2) Q! ), where td is the 
time during which 7r-pulses are applied [Fig. 2(b)]. It is 
difficult, however, to accurately determine a by directly 
fitting to this form. In contrast, we find that the scaling 
relation T2 ~ (n^) 1 very accurately describes the data 
irrespective of the value of a used to extract T2. The 
resulting 7 can then be related to a and other quanti- 



ties of interest within specific noise models. For CPMG 
with even n^, the scaling relation T2 oc (n^) le yields 
7 e = 0.72 ± 0.01, using T 2 values extracted using any 
a in the range 2 to 5. A model of dephasing due to 
a power-law spectrum of classical noise, S(u) ~ u; _/3 , 
leads to a scaling relation in the number of 7r-pulses, with 
the exponent of the power law, /?, related to the scaling 
exponent by the simple relation j3 = 7 e /(l — 7e)- For 
the present experiment, 7 e = 0.72 thus yields /3 = 2.6. 
Further support for a power-law form for S(uj) is found 
by comparing experimental and theoretical dependences 
T2(n 7r ) for CPMG with both even and odd rijc as well 
as CDD pulse sequences. This model also gives the sim- 
ple relation a = /3 + 1 connecting the noise spectrum 
and the decoherence envelope exponent. The resulting 
value, a = 3.6 ±0.1 is thus determined with considerably 
greater accuracy than can be obtained from direct fits to 
the coherence envelope data. 

The lateral double quantum dot investigated was de- 
fined by Ti/Au depletion gates patterned using electron 
beam lithography on a GaAs/Alo.3Gao.7As heterostruc- 
ture with a two dimensional electron gas (density 2 x 10 15 




FIG. 1: (Color online) (a) Micrograph of lithographically 
identical device with dot locations depicted. Gate voltages, 
Vr(l) 3 set the charge occupancy of right (left) dot as well 
as the detuning of the qubit. An rf-sensor quantum dot is 
indicated on the right, (b) Double dot charge state mapped 
onto dc conductance change, Ag, with lettered pulse sequence 
gate voltages. Detuning axis is orthogonal to the (0,2)-(l,l) 
charge degeneracy through points E, S, and M. 
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m -2 , mobility 20 m 2 /V s) 100 nm below the surface. 
Measurements were performed in a dilution refrigerator 
with an electron temperature T e ~150 mK. The double 
quantum dot is operated as a spin qubit by first depleting 
the quantum dots to the last two electrons, then manip- 
ulating the charge occupancy of the two dots with high 
bandwidth plunger gates Vl and Vr along a detuning 
axis e [Fig. 1(b)]. In this work, the charge occupancy 
was manipulated between states (0,2) and (1,1), where 
(Nl,Nr) represent the charge in the left and right dots. 
Charge occupancy was determined by the conductance 
change, Ag, through a proximal sensor quantum dot, 
which in turn modulated the reflection coefficient of the 
radio- frequency (rf) readout circuit [23) [24]. 

The logical spin qubit subspace is spanned by the sin- 
glet (S = (|tl) - \if))/y/2)) and the m = triplet 
(T = (\U) + \if))/y/2)) states of two electrons. The 
m = ±1 triplet states were split off by a 750 mT mag- 
netic field applied in the plane of the electron gas, per- 
pendicular to the dot connection axis. A (0,2) singlet was 
prepared at point P, off the detuning axis, through rapid 
relaxation to the ground state, then moved to the sepa- 
ration point S in (1,1). Uncorrelated Overhauser fields in 
the two dots create an evolving Zeeman gradient, AB Z1 
that drives transitions between S and Tq. Single-shot 
readout was performed by moving to point M, where S 
can tunnel to (0,2) while Tq remains in (1,1). The re- 
flectometer signal was integrated for 600 ns per shot, 
averaged over 10 4 shots, and compared to voltage val- 
ues corresponding to S and Tq outcomes [25], yielding 
-Ps'(td)? the probability of singlet return. 

Coherence lost due to (thermally driven) evolution of 
AB Z can be partially restored using a Hahn echo by puls- 
ing at time tjj/2 to point E, where the exchange splitting 
between S and To drives a 7r-rotation about the x axis, 
changing the sign of the acquired phase. Returning to S 
for an equal time r^/2 cancels the phase acquired due to 
the low- frequency (u<2/td) end of the spectrum of fluc- 
tuations of AB Z [TTJ [12] . Dynamical decoupling using a 
series of 7r-pulses allows efficient removal of more of the 
low-frequency end of the noise spectrum [TTJ [12]. The 
CPMG sequence [26 , for example, uses evenly spaced 
gate pulses from point S to point E with a half interval be- 
fore the first and after the last 7r-pulse [Figs.[2jb,c)]. Con- 
catenated dynamical decoupling [9] (CDD) uses nonuni- 
formly spaced pulses to point E, where the k-th order 
sequence is determined recursively from the lower order 
one, with an additional 7r-pulse in the center of odd or- 
ders [Fig. |2jd)]. 

Singlet return probabilities Ps(tb) were measured for 
CPMG sequences with n n = 1, 2, 3, 4, 8, 16, and 32. 
Fits to P s = 0.5 + V/2exp[-(r D /T 2 ) a ], with visibility 
V, T2, and a as a fit parameter, were equally good for a 
between 2 and 4, as seen in Fig. 3(a). For this reason, 
though S(uj) is related to a, these fits give little infor- 
mation about the spectrum of the environment. Figure 
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FIG. 2: (Color online) (a) Energy level diagram along the 
detuning axis e. The (0,2) singlet was prepared at ep, followed 
by separation to es- 7r-pulses are performed at 6e, allowing 
for subsequent rephasing at es- Single-shot readout at cm 
using a proximal sensor dot. The exchange energy, Je, that 
drives the 7r-pulses at ce is indicated with a dashed line, (b- 
d) Schematics of detunings during CPMG and CDD pulse 
sequences with detuning points on the vertical axis. 

3(b), showing T2 as a function of for an even number 
of CPMG pulses, shows two remarkable features. First, 
values of T2 do not depend on the value of a used in 
the fits to Ps'(td). Second, T2 shows a power-law scaling 
T 2 = T^in^) 1 whose power, 7, also does not depend on 
the value of a used in the fits. 

To model these observations, we consider Gaussian 
noise affecting the energy splitting of the qubit, which 
leads to the off-diagonal (in the basis of | ti) an d 
||t)) elements of the qubit density matrix decaying as 
exp[-x(r D )], where 

X(m)= r^SH^l, (1) 
Jo n ^ 

with F(uTn) being the filter function determined by the 
sequence of 7r-pulses driving the qubit. For CPMG se- 
quence F(z) < (z/2n 7V ) 4: for z < 2n^^ i.e. F(z) strongly 
suppresses the low-frequency noise, while for large z and 
riir the filter function can be approximated [12 by a pe- 
riodic train (with period z p = 2i\n^) of square peaks of 
height h~2n 2 K and width AzttZp/irh<^z p . 

We find that the value of 7 e and the presence of an 
even-odd effect (EOE) in n n (i.e. 7 e 7^ 7 Q ) act as dis- 
criminators for several classes of S(u): (i) The case of 
< 7e < 2/3 and absent EOE is compatible with a model 
of S(lj) ~ uo~P (over a range of uo roughly bounded by the 
minimal and maximal values of n^/ro), with < j3 < 2. 
In this case 7 e = f3/(l + j3) and a = j3 + 1. One exam- 
ple is the case of Ornstein-Uhlenbeck noise [28] (having 
Lorentzian S(u) with uo~ 2 tail typically dominating the 
decoherence under dynamical decoupling) , where 7 = 2/3 
was confirmed by experiments on the NV center [4 . (ii) 
When an EOE is present, 7 e = 2/3 suggests a hard cut- 
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FIG. 3: (Color online) (a) Experimental singlet return prob- 
abilities as a function of time for CPMG with — 2, 4, 8. 
Fits to P s (m) = 0.5 + y/2exp(-(r D /T 2 ) a ), with a con- 
strained to 2 (solid curves), 3 (dotted curves) and 4 (dashed 
curves) [27]. It is difficult to determine a from these fits, 
(b) Extracted T2 for even-n^ CPMG sequences for a con- 
strained to 1, 2, 3, 4, and 5. Circle size proportional to x 2 
goodness of fit of Ps(tv) in (a). A power- law fit to the form 
ln(T 2 even (a)) = ln(T 2 °) + 7 e ln(n w ), shown for a=3 (dashed 
line) gives 7 e = 0.72. The fit value 7 e depends only weakly 
on a in the range 2 — 5 (inset). The weighted average over 
a = 2-5 yields 7e = 0.72 ±0.01. 

off in S(lj) at u c < 2/T2, in which case j Q = 1. (m) For 
cj c > 2/T2, i.e. for larger cj c or larger (leading to longer 
T2), the EOE disappears and 7 tends to 1 [14 . (iv) Fi- 
nally, the presence of the EOE and 2/3 <7 e < 1 indicate 
S(u) ~ with /3>2. 

Experimentally, we find 7 e = 0.72 for even number of 
CPMG pulses, and = 1 not along the scaling line, in- 
dicating an EOE. We conclude that scenario (iv) applies, 
namely S(uj) = A^/ojt 3 with 2 < f3 < 3. Using Eq. (I) 
and the CPMG filter function gives in the large-n^ limit 

X (t) « (Ar D )^ (-1 + ^) , (2) 

with a^S2+/3/ / 7r 2 (27r) /3 , where ^s = J2T=i(^ ~ \)~ 5 ^ anc ^ 
for odd (even) n n we have b Q « [327r(3 — /3)] _1 (6 e « 
[128tt(5-/3)]- 1 « bo/10), i.e., the fc/ra* term is negligible 
for even n^, while it gives a significant correction for 
small, odd n^. The EOE comes from the difference in 
the \ow-z behavior of the CPMG filter functions, which 
for z < 1 behave as F(z) ~ z^j2 h n% (z 6 /2 7 nl) for odd 
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FIG. 4: (Color online) T 2 for all measured for CPMG, 
extracted using a = 2, 3, and 4 (circles). Circle size for each a 
proportional to \ 2 goodness of fit. Theory (black solid curve) 
for p — 7e/(l — 7e) = 2.6 and integration of Eq. |l]) with 
the CPMG filter functions. Note that Eq. ^ captures the 
even/odd effect quantitatively for small n^. 

(even) n n . For f3>2 this leads to different contributions 
of very low uo to the integral in Eq. 0. For even n n , we 
find that x( t d) approximately reduces to (Arj^) 1+ ^a/n^ 
from which we obtain the f3 < 2 result of j e = /3/(l + /?) 
in this case as well. 

Assuming this form of S(lj), fits to the even 
[Fig. 3(b)] yield /3 = 2.6 and A' 1 = 3.6 /is. Using these 
two parameters we calculate odd-n^ values for T2 by nu- 
merically integrating Eq. 0. As shown in Fig. |4j the 
obtained value of T2 is in good agreement with the mea- 
sured value for n n = 1 (Hahn echo). We note that the 
large scaling of T2 is due to the behavior of S(u) 
at cj>7rn 7r /T2, which is ~ 0.3(n 7r ) 0,28 /is -1 here, while 
the EOE at small is due to behavior at uo < I/T2, 
which is ~ 0.15 /is -1 . The consistency between small- 
and large-n^ data indicates that S(u) ~ u~ 2 - 6 over this 
range of frequencies (i.e., uj/2tt ~ 10 — 100 kHz). The 
EOE behavior at low can be fit within scenario (ii) 
using S(lj) = A 3 /uj 2 with A~ x ~l /is and u; c ~0.08 /is -1 . 
However, this scenario crosses over to (Hi) for > 5, 
where 7 tends to 1. The resulting large-n^ behavior, 
T2 ^riir x 7 /is, departs significantly from the > 8 data 
in Fig. 4. 

Using S(uj) = AP +1 /uu^ with parameters A and /? fixed 
from the even-n^ CPMG fit, we can calculate the ex- 
pected dependence of T2(n 7r ) for the CDD pulse sequence 
using the known filter functions [12]. For 77,^=5, 10, and 
21 we get good agreement between the calculated and 
measured T 2 [Fig. 5]. For 77^=42, the experimental T 2 is 
shorter than predicted by theory, possibly reflecting an 
accumulation of errors for such a large number of pulses. 
Note that CDD was shown to be robust to pulse errors 
[29] only in the case of two-axis control (i.e. when the 
7r-pulses are about x and y axes alternately) and for a 
quasi-static bath. 

Summarizing, the measurements of qubit decoherence 
under dynamical decoupling with the CPMG pulse se- 
quence have been used to reconstruct the crucial features 
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FIG. 5: (Color online) (a) Experimental singlet return 
probabilities for orders 1 through 6 CDD (symbols) along 
with fits to Ps(td) = 0.5 + V/2exp(-(TD/T 2 ) a )) using 
P = 7e /(l- 7e ) = 2.6 (solid curves), (b) Extracted T 2 
for using different values of a (circle size proportional to \ 2 
measure of goodness of fit). Theory (solid curve) based on 
P — 7e/(l — 7e) = 2.6 and integration of Eq. with appro- 
priate CDD filter function [12] . 



of the spectral density of noise dephasing the qubit. Us- 
ing the data for even of CPMG we have been able to 
estimate the magnitude of noise and its functional form. 
The reconstructed spectral density of noise allows us to 
calculate the expected decoherence signal for other pulse 
sequences, and this calculation agrees with the CDD se- 
quence measurements for as well as the odd-n^ CPMG 
data. We have shown that instead of fitting the exact 
functional form of the coherence decay function, an anal- 
ysis of the scaling of the measured T2 time with the num- 
ber of applied pulses allows for a clearer understanding 
of the system. We cannot say at this point whether the 
observed value j e = 0.72 is characteristic of Overhauser- 
dominated dephasing in general, or just our particular 
combination of noise sources. 
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